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Absorption of an electron by a dielectric wall 
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We introduce a method for calculating the probability with which a low-energy electron hitting 
the wall of a bounded plasma gets stuck in it and apply the method to a dielectric wall with positive 
electron affinity smaller than the bandgap using MgO as an example. In accordance with electron 
beam scattering data we obtain energy-dependent sticking probabilities significantly less than unity 
and question thereby for electrons the perfect absorber assumption used in plasma modeling. 

PACS numbers: 68.49.Jk, 79.20.Hx, 52.40.Hf 


The interaction of electrons with surfaces plays a key 
role in Mplied science. Various methods of surface anal¬ 
ysis [iRsI are based on it as well as a number of materials 
processing techniques In these applications the elec¬ 
tron energy is above 100 eV and backscattering and sec¬ 
ondary electron emission, the physical processes involved, 
are sufficiently well understood iEi. The situation is 
different for electron-surface interaction at energies below 
lOOeV, as it occurs in dielectric barrier discharges [il¬ 
ls I, dusty plasmas [H-ill , Hall thrusters [2l|, 22 , and 


electric probe measurements 23|. Much less is quanti¬ 
tatively known about it, especially at very low energies, 
below 10 eV. For instance, the probability with which 
a low-energy electron gets stuck in the wall after hit¬ 
ting it from the plasma is unknown. In the modeling of 
bounded plasmas [ii-il it is assumed to be close to 
unity, implying for electrons the wall is a perfect ab¬ 
sorber [2g|, irrespective of the material. Since electron 
absorption and extraction (by charge-transferring heavy 
particle collisions) control the wall potential, and hence 
the plasma sheath, which in turn affects the bulk plasma, 
the electron sticking probability is a crucial parameter. 
That its magnitude matters and should be known pre- 
cisely has been recognized most clearly by Mendis [20j in 
connection with grai n charging in dusty plasmas but the 
theoretical work [30l [31| he refers to is based on classical 
mechanics not applicable to electrons. 


In this work we apply quantum mechanics to calcu¬ 
late the probability with which a low-energy electron is 
absorbed by a surface. We couch the presentation in a 
particular application: The calculation of the electron 
sticking probability at room temperature for a dielectric 
wall [32j| with positive electron affinity x smaller than 
the bandgap Eg. But the approach is general and can be 
also applied to other cases. It utilizes two facts noticed 
by Cazaux [l[: (i) low-energy electrons do not see the 
strongly varying short-range potentials of the surface’s 
ion cores but a slowly varying surface potential and (ii) 
they penetrate deeply compared to the lattice constant 
into the surface. The scattering pushing the electron back 
to the plasma occurs thus in the bulk of the wall suggest¬ 
ing the probability for the electron to get absorbed by it 
to be the probability for transmission through the wall’s 


surface potential times the probability to stay inside the 
wall despite of internal backscattering. 

The variation of the potential across a floating dielec¬ 
tric wall can be calculated self-consistently [SSj. For 
X > 0 it gives rise to an energy barrier for electrons whose 
height on the plasma side is the Coulomb energy an 
electron has to overcome to reach the wall. On the solid 
side the height is x since for y > 0 electrons entering 
the wall and making up its charge are surplus electrons 
occupying the wall’s conduction band [33| . For the stick¬ 


ing probability it is the kinetic energy of the approaching 
electron in the vicinity of the wall which matters. It suf¬ 
fices therefore to model the wall by a three-dimensional 
potential step with height x ^'iid electron mass mismatch 
rrie = rn* jm^ < 1, where to* is the effective electron mass 
in the conduction band of the wall and rrie is the electron 
mass, as illustrated in Fig. [TJi. More refined treatments 
are possible. The shape of the step however is largely 
irrelevant for the energies we consider. For AI2O3, for in¬ 
stance, the attenuation length for an electron with a few 
eV is around 200 A and thus much larger than the 
lattice spacing setting the scale on which the step varies. 

Assuming as in Fig. [1^ the wall (plasma) to occupy 
the z < 0 (z > 0) half space, measuring energy in Ryd¬ 
bergs from the bottom of the conduction band, that is, 
setting Ech = Hw — X = 0, length in Bohr radii, and 
mass in electron masses, the transmission probability for 
a plasma electron hitting the wall, and having thus ki¬ 
netic energy A — x > 0, is given by [s^ 


T(E.f) = A-’"'*’’ 


{rriek +py 


( 1 ) 


with k = y/E — x^ and p = y/frieE p the z—components 
of the electron momenta outside and inside the wall. 
In IT]) the signs of k and p are always the same. We 
can thus define the direction cosines ^ and 77 referenced, 
respectively, to the electron momenta outside and in¬ 
side the wall, by their absolute values: ^ = |cos/3| and 
p = I cosd| (see Fig. [T]d for the definition of /3, 0, and 
E). This choice is dictated by the modeling of internal 
backscattering presented later. Since the potential varies 
only perpendicularly to the interface the lateral momen¬ 
tum K is conserved. Together with energy conservation. 
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E = x + k"^ + + ^^)/we, this yields 

E-X, 


1 - 7 ^^ = 


rripE 


(1-a 


( 2 ) 


connecting rj and From (21) follows that an electron 
from the plasma approaching the wall with kinetic energy 
E — X = > 0 enters it only when 


e>^c = 


for X < E < Eq 


/l-f^ for E>Eo 


( 3 ) 


with Eq = x/(l — me). For ^ the electron is in an 
evanescent wave with p'^ < 0, once it is inside the wall, 
and is thus totally reflected [Sg. This fact, unnoticed so 
far in connection with dielectric walls, has strong impli¬ 
cations for the sticking probability. In addition, Eq. (ED 
may force an electron entering the wall from the plasma 
to instantaneously acquire inside the wall a perpendicu¬ 
lar kinetic energy p^/me < X- For Trie < 1, applicable 
to MgO, the material we use for illustration because of 
the availab ility of data from electron beam scattering ex¬ 
periments 371, this happens when ^ < -^/l — ml. If such 
an electron cannot gain energy by inelastic scattering, as 
it is the case for dielectric walls at room temperature, it 
will have no chance to come back to the plasma since it 
will not be able to overcome the potential step x- 

The transmission probability captures only the ballis¬ 
tic aspects of electron absorption by the wall. Once the 
electron is inside the wall it is subject to elastic and in¬ 
elastic internal scattering, which may push the electron 
back to the interface and, after successfully traversing 
the surface potential in the reverse direction, back to the 
plasma (see Fig.[T^ for an illustration). To take this pos¬ 
sibility into account we define the sticking probability 
S{E,^) as the probability of an electron hitting the wall 
from the plasma with energy E and direction cosine ^ 
not to return to it after entering the wall and suffering 
backscattering inside it. Hence, 


S{E,0=r{E,0[l-£{E,0] 


( 4 ) 


with the conditional probability 

jlliJ^,dE'p{E')Q{Erj{0\EW)TiE',av')) 

£{E,Cj = -g- 

/o dfj' dE'p(E')Q(Ep(OlE'rj') 

( 5 ) 

for the electron to escape from the wall after at least 
one backscattering event. The lower integration limits in 
the numerator, r/rnin = \fxjE and EE^^ = ensure 

that only events are counted for which the perpendicular 
post-collision energy /me > Xi p{E) = I2{2ttY 

is the conduction band density of states, and Q{E'q\E'ri') 
is proportional to the probability for an electron, pene¬ 
trating the wall in a state with energy E and direction 



FIG. 1. (color online) (a) Potential energy step of height x 
modeling the wall and a scattering trajectory bringing an elec¬ 
tron entering the wall back to the plasma. Half circles are the 
moduli of the electron momenta outside and inside the wall 
at the entrance {E) and exit [E') point connected by phonon 
emission events. Also shown is the Coulomb barrier [/w the 
electron has to overcome to reach the wall, the energy dis¬ 
tribution f{E) it may have, and a realistic surface potential 
(red dashed curve). The tilt of the potential step due to the 
wall charge responsible for [/w is not shown, (b) Illustration 
of Eqs. (|4])--((5| defining S{E,^). The potential step leads to a 
transmission probability T while the emission of phonons to 
a conditional backscattering probability £ proportional to a 
quantity Q, obtained from the invariant embedding approach, 
that is, by adding an infinitesimally thin layer to the wall and 
requiring Q to be invariant against it [T^. Notice the defi¬ 
nition of the angles and energies characterizing the electron. 
Inside (outside) the wall the angles, measured with respect to 
the surface normal h, are denoted 9 {(3) and 9' [P'). 


cosine p, to backscatter after an arbitrary number of in¬ 
ternal scattering events towards the interface in a state 
with energy E' and direction cosine p' (see Fig. nt). Since 
Q{Ep\E'p') describes backscattering, 0 < 9 < 7r/2 and 
tt/2 < 9' < tt, implying 0 < p,p' < I as p' = |cos0'|. 
The energy integrals in (O anticipate that for a dielectric 
wall at room temperature the electron cannot gain energy 
and p{p) and its inverse P,(p) are defined by Eq. ([2). 

To calculate Q{Ep\E'p') we use the invariant em¬ 
bedding approach, originally applied to the electron 
backscattering problem by Dashen [l^ and revitalized 
recently by Vicanek Q and Glazov and coworkers d, . 
It is illustrated in Fig. [T]d. Summing up the four paths 
(l)-(4), requiring Q to be invariant against the change 
it induces, and taking into account that Q does not de¬ 
pend on the azimuths of the electron’s initial and final 
propagation direction measured with respect to the sur¬ 
face normal fi yields a nonlinear integral equation for 
Q{Ep\E'p') [13,115 ■ Since at low energy most scatter¬ 
ing processes in a solid are forwardly peaked [15 we lin¬ 
earize the equation by an expansion of Q{Ep\E'p') in the 
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number of the rare backscattering events. For dielectric 
surfaces at room temperature and 0 < x < Fig scat¬ 
tering at low energy arises mainly from the emission of 
optical phonons with energy lo. Since the electron can¬ 
not gain energy by scattering it can loose at most the 
energy it initially had when entering the wall. Expand¬ 
ing Q{Erj\E'rj') also in the number of forward scattering 
events yields thus a double series which terminates after 
a finite number of terms. From the golden rule transi¬ 
tion rate per time 0 and the material parameters of 


MgO 14^ follows that backward scattering due to 
emission of an optical phonon is two orders of magnitude 
less likely than forward scattering. Writing 


Q{Er^\E'i^') = Y.Y. vW)SiE - E'- <) (6) 

n—1 m—0 


with = (n+m)uj we can thus truncate the summation 
already after a single backward scattering event, that is, 
after n = 1 leaving Mtot = [E /wj — 1 forward scattering 
events at most. Introducing kernels 

K^{Erj\E'r,') = [{E + E'T 

-4EE'{1-7j^){1-i^'^)]-^^\ (7) 


for forward (-I-) and backward (—) scattering and a func¬ 
tion n(E) = arcosh( ^EjiE) jE, all three arise in the 
course of the calculation from the transition rate per 
length, the input driving the integral equation for Q 


most directly obtained by dividing the transition rate per 
time by the electron velocity prior the emission event, the 
expansion coefficients satisfy 


Qln{E\r]\r]') = F^{E]r]\r]')Ql^_^{E - u:]r]\r]') 

+ Qln-i{E]'n\'q')Grn{E-mw\'q\vi') ( 8 ) 


with TO = 1 ,..., Mtot, and 


F (F-r,\E\ K+{E\E-u:]ri)ip{E-u}) 

ri'\l{E)+r]\l{E-{m+l)uj) ' 
^ VP{E)K+{E\E 

Gm\E,Tj\r] ) — ml TP I ^ I TllJP \ ' 

r]'ll(E + moj) + r]ll{E — uj) 
The initialization of the recursion ([5]) is given by 


QliE;vW) = ^ 


rj'K {Eri\E — ujrj') 


r]'Ii{E) -\- r]Ii{E — uj) 


(9) 

( 10 ) 

( 11 ) 


Forward scattering is encoded into the function 
K+{E\E'-i^) = ^ldfiK+{Er]\E’fi) = dfiK+{Ef]\E’ri), 
where the second identity utilizes the symmetry of 
K'^{Er]\E'rj') with respect to interchanging 77 and 77 '. It 
can be obtained analytically and enters the formalism be¬ 
cause K^{Eri\E''q') is strongly peaked for 77 = 77 ', as can 
be seen from Eq. 0 , and noticing that K'^{Er]\E'r]') en¬ 
ters always with E' = E — ui. Integrals over the direction 



Energy in eV Energy in eV 

FIG. 2. (color online) Electron sticking probability for MgO 
obtained from Eq. On the left is shown S{E,^) for all di¬ 
rection cosines ^ and energies E up to the bandgap Eg. Total 
reflection occurs in the white region. Below the yellow dotted 
line inelastic backscattering has no effect on S{E,^). On the 
right S{E,^) (solid curves) and T{E,^) (dashed curves) are 
plotted as a function of E for representative The material 
parameters are Trie = 0.4, x = leV, Eg = 7.8 eV [^, and 
LO = 0.1 eV [ 4 ^. The value for me is strictly applicable only 
at the band edge. Away from it me depends on momentum. 
From the bandstructure [dj we expect our results reliable up 
to 5eV enough for a comparison with experiment (see Fig. El. 


cosine containing K^{Eri\E'r]') can thus be handled by 
a saddle-point approximation. Physically this means the 
directional change due to forward scattering is negligi¬ 
ble. Only backward scattering changes the direction (see 
Fig. [Dr). We thus end up with a model similar in spirit to 
the Oswald, Kasper and Gaukler model ^ for backscat¬ 
tering of electrons from surfaces due to multiple elastic 
scattering. The strong forward scattering comes from the 
matrix element for the interaction of an electron with an 
optical phonon which is at low electron density inversely 
proportional to the momentum transfer squared (40l |. 

Substituting expansion ([ 6 ]) for Q{Ef]\E'ri') in dS]) sets 
E' to El^ = E — and replaces the energy integrals 
by sums over to running in the numerator up to to = 
Mopen = [(.Erf — x)/(rf^)\ — 1 and in the denumerator 
up to TO = Mtot- Inserted into Eq. (jj]) it finally gives the 
electron sticking probability S{E, at room temperature 
for a dielectric wall with 0 < x < Eg and a clean, lateral 
momentum conserving interface. 

In Fig. ID we present the results for MgO obtained from 
Eq. (ID). Material parameters are given in the caption. 
The left panel shows the angle-resolved sticking probabil¬ 
ity S{E,^) over the whole range of direction cosines and 
energies up to the bandgap, above which electron-hole 
pair generation across the bandgap has to be taken into 
account as an additional, Coulomb-driven backscattering 
process. It can be treated in the same spirit, in particular 
the truncation after one backscattering event will still be 
possible, but the recursion will then contain energy inte¬ 
grals. The white area in the plot for S{E, ^) indicates the 
region in the [E, ^)-plane where total reflection occurs. It 
is large because nie is for MgO significantly smaller than 
unity. Below the dotted yellow line, ^ = 1 /I — TOe, in¬ 
elastic backscattering due to emission of optical phonons 
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is irrelevant for sticking because the perpendicular en¬ 
ergy of the electron drops, upon entering the wall, below 
the potential step x due to conservation of total energy 
and lateral momentum. It is hence already confined by 
transmission. Only above the dashed yellow line inelastic 
backscattering has a chance to bring the electron, once it 
is inside the wall, back to the interface and, after a suc¬ 
cessful reverse transmission through the surface poten¬ 
tial, back to the plasma. As a result, S{E,^) = T{E,^) 
for ^ < y/1 — rrie and S{E, < T{E, for ^ — rue- 

This can be more clearly seen in the right panel, where 
S{E,^) (solid curves) and T{E,^) (dashed curves) are 
plotted as a function of E for representative 

The sticking probabilities in Fig. [5] are for ^ < Cc 
strongly affected by total reflection caused by the mass 
mismatch and conservation of energy and lateral momen¬ 
tum. For the latter in-plane homogeneity of the interface 
is crucial. In reality imperfections may destroy it and 
lateral momentum may thus not be conserved. To take 
this possibility into account we include hard core poten¬ 
tials to mimic interfacial scattering centers and adopt 
the modeling of ballistic electron emission spectroscopy 
by Smith and coworkers [i^ to the calculation of S. For 
a disordered interface we then obtain 

+ I' ’ ( 12 ) 

where £{E,^) is given by ([5]) with T{E,^) replaced by 

The parameter (7 is a fit parameter proportional to the 
density of scattering centers. It controls the elastic scat¬ 
tering at the interface. In the weak scattering limit, 
C —>■ 0 and S' —>■ S'. We thus recover the sticking prob¬ 
ability for the clean interface. In the strong scattering 
limit, C —>■ oo and we get the sticking probability for the 
totally disordered interface for which C is irrelevant. 

Results for MgO obtained from Eq. m are shown 
in Fig. [3] for ^ = 1 (normal incident) and C = 0,1,2, 
and 00 . We also plot experimental data from electron 
beam scattering [37|. The agreement between calculated 
and measured probabilities is astonishing indicating our 
approach captures essential aspects of electron absorp¬ 
tion by the wall at low energies. Dashed curves show 
T{E, 1), the sticking probability in the absence of inter¬ 
nal backscattering. For C = 0, S{E, 1) deviates strongly 
from T(E, 1) (black curves), whereas for C = oo the two 
quantities approach each other (blue curves). The reason 
is the angle-averaging at the disordered interface which 
lessens, for fixed the impact of internal backscattering 
compared to the knock-out of ^-values by total reflec¬ 
tion signalled by the kink at E = Eq. How close the 



FIG. 3. (color online) Electron sticking probability S{E-, 1) 
for normal incident onto a MgO wall obtained from Eq. (HD 
(solid curves). We show data for C = 0 (black), 1 (red), 2 
(green), and oo (blue) with T{E\ 1) also included for C = 0 
and (7 = 00 (dashed curves). Symbols are data from electron 
beam scattering experiments, where two methods, denoted 
potential adjusted method (pam) and potential subtraction 
method (psm), have been employed to determine the energy 
of the electron [^. Material parameters as in Fig. [3] 

two curves come to each other, that is, how ineffective 
internal backscattering becomes depends on me and y. 

To summarize we presented a quantum-mechanical 
method to calculate the sticking probability S for a low- 
energy electron hitting the wall of a bounded plasma. 
Using the electron’s large penetration depth the method 
expresses S as the probability E for transmission through 
the wall’s surface potential times the probability 1 — 8 
to stay inside it despite of internal backscattering. For 
MgO we got in agreement with electron beam scattering 
data S' < 1, indicating S ~ 1, employed in plasma 
modeling, may not always be justified. Using S{EE) 
from our approach instead of S ~ 1, for instance, in the 
theory of grain charging in dusty plasmas would reduce 
the charge of micron-sized MgO particles by more than 
ten percent. The method is flexible and numerically effi¬ 
cient. Various scattering processes, including interfacial 
scattering due to disordered interfaces, can be treated. 
It could thus guide an experimental effort to develop 
a realistic description of electron-wall interaction as it 
occurs in many applications of bounded plasmas. 

We acknowledge support by the Deutsche Forschungs- 
gemeinschaft through the Transregional Collaborative 
Research Center SFB/TRR24. 


[1] J. Cazaux, J. Appl. Phys. Ill, 064903 (2012). 

[2] F. Salvat-Pujol and W. S. M. Werner, Surf. Interface 
Anal. 45, 873 (2013). 

[3] W. S. M. Werner, Surf. Interface Anal. 31 , 141 (2001). 

[4] A. Klassen, A. Bauereiss, and C. Koerner, J. Phys. D: 
Appl. Phys. 47 , 065307 (2014). 















5 


[5] L. G. Glazov and S. Tougaard, Phys. Rev. B 68, 155409 
(2003). 

[6] E. Schreiber and H.-J. Fitting, J. Electron Spectros. Re- 
lat. Phenom. 124, 25 (2002). 

[7] A. Dubus, A. Jablonski, and S. Tougaard, Progr. Surface 
Science 63, 135 (2000). 

[8] M. Vicanek, Surface Science 440, 1 (1999). 

[9] S. L. Dudarev, P. Rez, and M. J. Whelan, Phys. Rev. B 
51, 3397 (1995). 

[10] A. L. Tofterup, Phys. Rev. B 32, 2808 (1985). 

[11] K. Kanaya and S. Okayama, J. Phys. D: Appl. Phys. 5, 
43 (1972). 

[12] R. Dashen, Phys. Rev. 134, A1025 (1964). 

[13] M. Bogaczyk, R. Wild, L. Stollenwerk, and H.-E. Wag¬ 
ner, J. Phys. D: Appl. Phys. 45, 465202 (2012). 

[14] R. Tschiersch, M. Bogaczyk, and H.-E. Wagner, J. Phys. 
D: Appl. Phys. 47, 365204 (2014). 

[15] F. J. J. Peelers and M. G. M. van de Sanden, Plasma 
Sources Sci. Technol. 24, 015016 (2015). 

[16] P. Tolias, Plasma Phys. Control. Fusion 56, 123002 
(2014). 

[17] I. Richterova, M. Beranek, J. Pavlu, Z. Nemecek, and 
J. Safrankova, Phys. Rev. B 81, 075406 (2010). 

[18] O. Ishihara, J. Phys. D: Appl. Phys. 40, R121 (2007). 

[19] V. E. Fortov, A. V. Ivlev, S. A. Khrapak, A. G. Khrapak, 
and G. E. Morfill, Phys. Rep. 421, 1 (2005). 

[20] D. A. Mendis, Plasma Sources Sci. Technol. 11 , A219 
( 2002 ). 

[21] A. Dunaevsky, Y. Raitses, and N. J. Fisch, Phys. Plasmas 
10, 2574 (2003). 

[22] S. Barral, K. Makowski, Z. Peradzyhski, N. Gascon, and 
M. Dudeck, Phys. Plasmas 10, 4137 (2003). 

[23] M. Lapke, T. Mussenbrock, R. P. Brinkmann, C. Schar- 
witz, M. Boeke, and J. Winter, Appl. Phys. Lett. 90, 
121502 (2007). 

[24] M. Shihab, A. T. Elgendy, I. Korolov, A. Derzsi, 
J. Schulze, D. Eremin, T. Mussenbrock, Z. Donko, and 
R. P. Brinkmann, Plasma Sources Sci. Technol. 22, 
055013 (2013). 

[25] L. Stollenwerk, S. Amiranashvili, J.-P. Boeuf, and H.-G. 
Purwins, Eur. Phys. J. D 44, 133 (2007). 


[26] H. Kersten, H. Deutsch, and G. M. W. Kroesen, Int. J. 
Mass Spectr. 233, 51 (2004). 

[27] R. Brandenburg, V. A. Maiorov, Y. B. Golubovskii, H.- 
E. Wagner, J. Behnke, and J. F. Behnke, J. Phys. D: 
Appl. Phys. 38, 2187 (2005). 

[28] D. Uhrlandt, M. Schmidt, J. F. Behnke, and T. Binde- 
mann, J. Phys. D: Appl. Phys. 33, 2475 (2000). 

[29] Y. L. Al’pert, A. V. Gurevich, and L. P. Pitaevskii, 
Space physics with artificial satellites (Consultants Bu¬ 
reau, New York, 1965). 

[30] T. Umebayashi and T. Nakano, Publ. Astron. Soc. Japan 
32, 405 (1980). 

[31] D. Hollenbach and E. E. Salpeter, J. Chem. Phys. 53, 79 
(1970). 

[32] F. Bechstedt, Principles of surface physics (Springer Ver- 
lag, Berlin, 2003). 

[33] R. L. Heinisch, F. X. Bronold, and H. Fehske, Phys. Rev. 
B 85, 075323 (2012). 

[34] T. W. Hickmott, J. Appl. Phys. 36, 1885 (1965). 

[35] C. M. Wu and E. S. Yang, Solid-State Electronics 22, 
241 (1979). 

[36] T. K. Gaylord and K. F. Brennan, J. Appl. Phys. 65, 
814 (1989). 

[37] C. J. Cook and W. J. Fredericks, J. Chem. Phys 36, 608 
(1962). 

[38] L. G. Glazov and I. Pazsit, Nucl. Instr. and Meth. B 256, 
638 (2007). 

[39] B. K. Ridley, Quantum processes in semiconductors 
(Clarendon Press, Oxford, 1999). 

[40] B. K. Ridley, J. Appl. Phys. 84, 4020 (1998). 

[41] Y.-S. Kim, S.-H. Yoon, S.-G. Ahn, C.-R. Hong, and 
H. Yang, Electron. Mater. Lett. 4, 113 (2008). 

[42] A. R. Oganov, M. J. Gillan, and G. D. Price, J. Chem. 
Phys. 118, 10174 (2003). 

[43] R. Oswald, E. Kasper, and K. H. Gaukler, J. Electron 
Spectros. Relat. Phenom. 61, 251 (1993). 

[44] Y.-N. Xu and W. Y. Ching, Phys. Rev. B 43, 4461 (1991). 

[45] D. L. Smith, E. Y. Lee, and V. Narayanamurti, Phys. 
Rev. Lett. 80, 2433 (1998). 



